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A comprehensive phenomenological analysis of a two Higgs doublet model, with flavor-changing 
scalar currents at the tree-level, called model 111, is presented. Constraints from existing experi- 
mental information especially on AF — 2 processes are systematically incorporated. Constraints 
\^ I emerging from rare B-decays, Z hb, and the p-parameter are also examined. Experimental 

^\ . implications for e"'"e~(^"'"/i~) tc + tc, t c^{Z,g), D'^-D'^, and B^-B^ oscillations, and for 

^\ ' e'^e~{Z) — > ba -\-bs are investigated and experimental effort towards these is stressed. We also 

emphasize the importance of clarifying the experimental issues pertaining to Z bb. 
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I. INTRODUCTION 



Q^ , Recently we have been examining various issues in a class of two Higgs doublet models (2HDM's) which allow 
C ' flavor changing neutral currents (FCNC's) at the tree level In this work we want to present a comprehensive 

analysis which gives the details of the analytical calculations we used and summarize the status of our knowledge 
regarding this type of 2HDM's. In particular we will examine the important constraints and derive quantitative 
bounds on the mass parameters and flavor changing (FC) couplings of the new scalar fields based on existing low 
. energy experiments. 

' FCNC's are naturally suppressed in the standard model (SM) because they are forbidden at the tree level. However, 
, at the one loop level, the amplitudes for the Feynman diagrams which generate FC processes tend to increase with 
II virtual quark masses. Due to the large disparity in masses of the up-type quarks, the GIM |^ suppression gets very 
Oh! effectively removed in loops of these (virtual) quarks. Therefore, FC transitions involving a pair of down-type quarks 
' get enhanced, due primarily to the presence of a top-quark in the loop, while FC transitions which involve up-type 
quarks are usually very small. This motivates the interest in processes like 6 — > s, ^7 instead of i ^ C7, Z or for 
j^a _ j^Q ^j^j j^Q _ j^o jj^ixing instead of — mixing. 

There are many ways in which the extensions of the SM lead to FC couplings at the tree level. For instance, as 
soon as we go from a theory with one doublet of scalar fields to a theory with two doublets, FCNC's are generated 
in the scalar sector of the theory. In this case, couplings such as tcj or tcZ can get enhanced too, in much the same 
manner as there is an enhancement of bsj in the SM. The interest in this class of implications is obvious: we could 
have a clear signal of new physics, since the SM prediction for any process involving a tcy or a tcZ vertex is extremely 
small. 

As is well known, a model with tree level FCNC's will also have many important repercussions for the AF — 2 
mixing processes such as if" — K^, B'^ — B^, and D° — l)". Indeed the measured size of the mixing amplitude for 
— K^^ known for a very long time, is so small that it places severe restrictions on the FC sector of extended models. 
This led Glashow and Weinberg ||] to propose an ad hoc discrete symmetry whose sole purpose was to forbid tree- 
level FCNC's to appear in models with more than one Higgs doublet. In particular, for the simplest case of 2HDM's, 
depending on whether the up-type and down-type quarks couple to the same or to two different scalar doublets leads 
to two versions of such models, called model I and model II. Both of them share with the SM the distinctive feature 
that they do not allow tree-level FCNC's. These models have been extensively studied and a multitude of interesting 
implications have been discussed in the literature. 

Our starting point for investigating tree-level FCNC's is primarily based on the realization that the top quark may 
be quite different from the lighter quarks. It may well be that the theoretical prejudice of the non existence of tree 
level FCNC's, based on experiments involving the lighter quarks, is not relevant to the top quark. This leads one 
to formulate a model which allows the possibility of large tree level FCNC's involving the top quark while it keeps 
the FCNC's of the lighter quarks, especially those involving the quarks of the first family, at a negligible level. A 
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rather natural way of implementing this notion is by taking the new scalar FC vertices to be proportional to the 
masses of the participating quarks at the vertex [^|-|5| . A hierarchy is then automatically introduced which enhances 
the top quark couplings while keeping the couplings of the lighter quarks at some order of magnitude smaller. The 
compatibility of any such assumption with the many constraints mentioned above clearly needs to be checked and the 
allowed region of the parameter space determined, in particular, scalar masses and couplings. This type of 2HDM is 
now referred to as model III ||lot] . 

Our first interest in model III was motivated by the idea of looking for top-charm production at an e^e~ and/or 
Ijf^ jjT collider We were intrigued by the possibility of a clear signal for the reaction e+e^ tc for which the 

SM prediction is extremely small. The reaction has a distinctive kinematical signature in a very clean environment. 
These characteristics, which are unique to the lepton colliders, should compensate for the lower statistics one expects 
compared to those at hadron colliders. In this paper we present the details of our calculations for this important 
process. 

After a brief overview of model III in Sec. 0, we will discuss e+e ^ tc + tc in Sec. CIJ and present some of the 
relevant formulae in detail in Appendix In Sec. |^ we will consider the rare decays t cj, Z and g; herein we will 
also compare our results with those existing in the literature for t cy and cZ in model III. Top-charm production 
at a fJ-'^ fJ.~ collider is briefly discussed in Sec. 0. In Sec. we consider the reper cuss ions of the AF = 2 mixing 
processes and extract the constraints that emerge on parameters of model III. Sec. VII discusses the i mpact of the 
experimental results for Br{B — > Xg^), the p-parameter and Z bb. We then examine in Sec. VIIl the physical 
consequences of the constrained physical model for the — B^ oscillations, the flavor-changing Z decay Z ^ bs + bs, 



and also some rare i?-decays. Sec. IX offers the outlook and the conclusions. 



II. THE MODEL 



A mild extension of the SM with one additional scalar SU(2) doublet opens up the possibility of flavor changing 
scalar currents (FCSC's) at the tree level. In fact, when the up-type quarks and the down-type quarks are allowed 
simultaneously to couple to more than one scalar doublet, the diagonalization of the up- type and down- type mass 
matrices does not automatically ensure the diagonalization of the couplings with each single scalar doublet. For this 
reason, the 2HDM scalar potential and Yukawa Lagrangian are usually constrained by an ad hoc discrete symmetry 
1^, whose only role is to protect the model from FCSC's at the tree level. Let us consider a Yukawa Lagrangian of 
the form 
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where 0i, for i = 1, 2, are the two scalar doublets of a 2HDM, while rj^^'^ and Cij'^ ^'^^ the non-diagonal matrices of 
the Yukawa couplings. Imposing the following ad hoc discrete symmetry 
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one obtains the so called model I and model II, depending on whether the up-type and down-type quarks are coupled 
to the same or to two different scalar doublets respectively | pT[ . 

In contrast we will consider the case in which no discrete symmetry is imposed and both up-type and down-type 
quarks then have FC couplings. For this type of 2IIDM, which we will call model III, the Yukawa Lagrangian for the 
quark fields is as in Eq. (0) and no term can be dropped a priori^ see also refs. WJlfl- 

For convenience we can choose to express „.,■,^-,.u^,. t,„„:„ „,,„i, j-i,„j- — i., 4-u — u,d 

the fermion masses, i.e. such that 
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The two doublets are in this case of the form 
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The scalar Lagrangian in the {H^, , H'^ , H^) basis is such that |ll2| , ^ : 
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1. the doublet (pi corresponds to the scalar doublet of the SM and to the SM Higgs field (same couplings and 
no interactions with and H^); 

2. all the new scalar fields belong to the 02 doublet; 

3. both and do not have couplings to the gauge bosons of the form H^'^ZZ or H^'^W^W~ . 

However, while is also the charged scalar mass eigenstate, , H^) are not the neutral mass eigenstates. 

Let us denote by (^°, hP) and the two scalar plus one pseudoscalar neutral mass eigenstates. They are obtained 
from (i/o, H^, H^) as follows 

= [{H° - v)cosa + sin a] 
h" =[-{H"-v)sma + H^cosa\ (5) 

where a is a mixing angle, such that for a — 0, (iJ", H^, H^) coincide with the mass eigenstates. We find it more 
convenient to express and as functions of the mass eigenstates, i.e. 

ij" = (i7"cosa- /i"sina) +v 

= {h° cos a + H° sin a) (6) 
i/2 ^ A" . 

In this way we may take advantage of the mentioned properties (1), (2), and (3), as far as the calculation of the 
contribution from new physics goes. In particular, only the 0i doublet and the r^Yj and r/^ couplings are involved in 
the generation of the fermion masses, while (j)2 is responsible for the new couplings. 

After the rotation that diagonalizes the mass matrix of the quark fields, the FC part of the Yukawa Lagrangian 
looks like 

^y^FC = lYA,LhUj,R + £.^JQ^,LhDJ,B. + h.c. (7) 

where (5i,L, Uj^r, and Dj r denote now the quark mass eigenstates and ifj'^ are the rotated couplings, in general not 

diagonal. If we define Vj^'l^ to be the rotation matrices acting on the up- and down-type quarks, with left or right 
chirality respectively, then the neutral FC couplings will be 

4noutral = [^L ) '^R ■ i^) 

On the other hand for the charged FC couplings we will have 

■ V 

Schargcd 'sncutral •'Ckm 

^charged ~ ^ckm ' '^neutral (9) 

where Vckm denotes the Cabibbo-Kobayashi-Maskawa matrix. To the extent that the definition of the ^^'^ couplings 
is arbitrary, we can take the rotated couplings as the original ones. Thus, we will denote by ^^j^ the new rotated 
couplings in Eq. (||), such that the charged coupHngs in Eq. (^) look hke ■ Fckm and V^ckm • ■ This form of 
the charged couplings is indeed peculiar to model III: they appear as a linear combination of neutral FC couplings 
multiplied by some CKM matrix elements. This is an important distinction between model III on the one hand, and 
models I and II on the other. As we will see in the phenomenological analysis this can have important repercussions 
for many different physical quantities. 

In order to apply to specific processes we have to make some definite ansatz on the (^-^ couplings. Many different 
suggestions can be found in the literature ^[0]. In addition to symmetry arguments, there are also arguments 
based on the widespread perception that these new FC couplings are likely to mainly affect the physics of the third 
generation of quarks only, in order to be consistent with the constraints coming from K^—K^ and B^—B'^. A natural 
hierarchy among the different quarks is provided by their mass parameters, and that has led to the assumption that 
the new FC couplings are proportional to the mass of the quarks involved in the coupling. Most of these proposals 
are well described by the following equation 
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FIG. 1. Example of one-loop contributions to the Ztc and 7tc vertices in model III. 



which basically coincides with what was proposed by Cheng and Sher Q. In this ansatz the residual degree of 
arbitrariness of the FC couplings is expressed through the Ay parameters, which need to be constrained by the 
available phenomenology. In particular we will sec how — K^ and B'^ — B^ mixings (and to a less extent D^ — D^ 
mixing) put severe constraints on the FC couplings involving the first family of quarks. Additional constraints are 
given by the combined analysis of the Br{B Xg'y), the p-parameter, and Rb, the ratio of the Z ^ bb rate to the 
Z-hadronic rate. We will analyze all these constraints in the following sections and discuss the resulting configuration 



of model III at the end of Sec. VIl 



III. TOP-CHARM PRODUCTION AT e+e" COLLIDERS 



The presence of FC couplings in the Yukawa Lagrangian of Eq. affects the top-charm production at both 
hadron and lepton colliders. In particular we want to study top-charm production at lepton colliders, because, as we 
have emphasized before [^]J^], in this environment the top-charm production has a particularly clean and distinctive 
signature. In principle, the production of top-charm pairs arises both at the tree level, via the s-channel exchange 
of a scalar field with FC couplings, and at the one loop level, via corrections to the Ztc and jtc vertices. We will 
consider in this section the case of an e+e" collider and in Sec. ^ that of a fi~^^J,~ collider. 

The s-channel top-charm production is one of the new interesting possibilities offered by a collider in studying 

the physics of standard and non-standard scalar fields (see Sec. ^ and refs. therein). However, it is not relevant for an 
e+e" collider, because the coupling of the scalar fields to the electron is likely to be very suppressed (see Eq. (p^). 
Therefore we will consider top-charm production via 7 and Z boson exchange, i.e. the process e+e~ — > 7*, Z* tc+ct, 
where the effective one loop 7tc or Ztc vertices are induced by scalars with FC couphngs (e.g. model III) [p^ . 

Let us write the one loop effective vertices Ztc and jtc in the following way 



*^ 16^ 



( A(^)7'' + 5(^)7^7^ + ^C<''^a^^-±- + w^""^ af^-" ^j^] tV, (11) 
V nit mt J 



where V = ^, Z and A^'^'^\ B^'^'^\ . . . denote the form factors generated by one loop corrections. A sample of the 
corrections one has to compute in model III isgiven in Fig. 1, for both the neutral and the charged scalar fields. The 
analytical expressions are given in Appendix y| where the details of the form factor calculation are also explained. 

In terms of the A^'^'^\ B^'''^\. . . form factors we can compute the cross section for e+e" — > 7*, Z* ^ ct + tc. The 
total cross section will be the sum of three terms 

ctot = + (Tz + <y-iZ (12) 
corresponding to the pure photon, pure Z and photon-Z interference contributions respectively. Thus 
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(1-/32) 

Re A'^Im A^ - Im ^''Re A^ + Re B''Im - Im B'^Re B^ 
(Re C'lm - Im C^Re + Re I^'^Im - Im ZJ'^Re 



(15) 



where Ue denotes the QED fine structure constant, Siv — s'md„, c„ = cos6'iv, s ~ is the center of mass energy 
squared, D{s) denotes the Z boson propagator and we have introduced 



/3^ = 1-^ 
s 



(16) 



We will consider the total cross section normalized to the cross section for producing pairs via one photon 

exchange, i.e. 



a{e^e tc + tc) 



cr(e+e 



(17) 



and normalized to ~ A = 1 (see Eq. (|1C|)). For the moment, we want to simplify our discussion by taking the 
same A for both ^^t ^^'^ ^ct- Moreover, we want to factor out this parameter, because it summarizes the degree of 
arbitrariness we have on these new couplings and it will be useful for further dis cuss ion. We will elaborate more about 
the possibility of considering different alternatives for the FC couplings in Sec. VII , after we present a comprehensive 
analysis of the constraints. 
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FIG. 2. 7?'7A* vs. ^ for case 1 (solid), case 2 (dashed), and case 3 (dot-dashed), with m* = 180 GeV. See Eq. (Q. 

As already discussed in Ref. |l[], we take mt — 180 GeV and vary the masses of the scalar and pseudoscalar fields 
in a range between 200 GeV and 1 TeV. Larger values of the scalar masses are excluded by the requirement of a 
weak-coupled scalar sector. The phase a does not play a relevant role and, as we discuss in Appendix H, in our 



5 



Decay SM Model I Model II Model III 



t->c7 ~5-10~^2 10-1^-10"" 10-"-10-^ 10-^^-10"^ 

t^cZ ~ 10"" 10""-10"" 10""-10"" 10"*-10"'' 

t^cg ~5-10"" 10""-10-'' 10""-10"* lO^^-lO"" 



TABLE I. Values of Br{t ^ C7), Br{t cZ), and Br{t eg) for mt ~ 180 GeV, in the SM and in the 2HDM's denoted 
as Model I, Model II, and Model III. Each range is obtained by varying rric, ruh, tua, tan/?, . . . over a broad region of the 
parameter space of the corresponding model, as explained in the text. For Model III, we have fixed \ij ~ A = 1 in the FC 
couplings. 



qualitative analysis we will set a = 0. In this case, as we can read in Eqs. (^) or (^, = + = h^, and 

= ^0 g^j^j ^Yie only new contributions come from and ^4°. In Fig. || we plot ff'^/X'^ as a function of ^/s for a 
sample of relevant cases, in which one of the scalar particles is taken to be light (Af; ~ 200 GeV) compared to the 
other two {Mh~l TeV), i.e. 
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where ruh and uia are the neutral scalar and pseudoscalar masses and rUc is the charged scalar mass respectively. 
We find that even with different choices of nih, mA and rric it is difficult to push i?*'^/A^ much higher than 10~^. 
Therefore the three cases illustrated in Fig. || appear to be a good sample to illustrate the type of predictions we can 
obtain for the rate for top-charm production in model III. 

From Fig. ||, we also see that going to energies much larger than ~ 400-500 GeV (i.e. ~ 2Mi) does not gain much 
in the rate and in this case R*'^/X'^ can be as much as 10~^. Since it is reasonable to expect 10''-10^ M^M^ events 
in a year of running for the next generation of e+e" colliders (/ £ ~ 5 x lO'^'^ cm~^sec~^) at ^/s = 500 GeV, this 
signal could be at the detectable level only for not too small values of the arbitrary parameter A. Thus we can expect 
experiments to be able to constrain A<^1, for scalar masses of a few hundred GeVs. 



IV. RARE TOP DECAYS: t C7, Z, g 

Starting from the form factors defined in Eq. ( pi] ) and given in Appendix |^, we can also easily derive the rates for 
rare top decays like t 07, t cZ and t eg. The study of rare top decays has been often emphasized in the 
literature |]l5| , p^ , p^| J^] , in particular, as a potential source of evidence for new physics. Indeed, as we can read from 
Table ^, these decays are extremely suppressed in the SM and they are quite small even in the 2IIDM's without tree 
level FCNC's (i.e. both in model I and in model II) This is due to a strong GIM suppression from the small 

value of the internal quark masses md^s,b as well as the large tree level rate for t bW. On the other hand, these rare 
top decays normally get enhanced in models with FGNC's and this motivates us to estimate their branching ratio 
in model III. From the experimental point of view the prospects for the three modes, t — *■ C7, cZ, and eg are quite 
different. In particular, t eg could be quite problematic for a hadron collider and the backgrounds will have to be 
considered before one can ensure that they do not represent a serious limitation. On the other hand, for the e~^e~ 
case background issues are less likely to be a serious problem even for t eg. 

In model III the modes t ^ cy and i — > have been previously considered and their rates are given by 

r(^->^7)=(^^(|c^l^ + m^) (19) 
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The rate for t eg can be written in an analogous manner, as, 

1 1 



where Cp = [N^ - l)/2N. 
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FIG. 3. Branching fraction for t 
mass Ms, when mt = 180 GeV. 
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C7 (solid), t cZ (dot-dashed), and t ^ eg (dashed) as a function of a common scalar 



The branching ratios reported in Table || are obtained by normalizing T{t 
simplicity, just to the main decay t bW rate, i.e. 



C7), T{t cZ) and T{t — > eg), for 
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The branching ratios for model I and model II are deduced from the analysis of Rcf. ||15|,|1^, with mt~180 GeV. The 
results for model III are obtained by varying the neutral and the charged scalar masses between 200 GeV and 1 TeV, 
assuming different patterns as explained in the preceding section (see for instance Eq. (p^)). In particular the upper 
bounds on the different branching ratios given in Table | are obtained by taking the scalars to have a common and 
relatively small mass. We also notice that in model III the results show a significant dependence on mt, such that the 
numbers in Table ^ change on the average by as much as an order of magnitude when mt is varied between 150 GeV 
and 200 GeV. This sensitivity to mt may become relevant when the experiments ever get to the point of being able 
to measure this type of rare top decay. Finally, in the FC couplings of Eq. (^^ we also take all the Ay parameters to 
be equal to A. In particular, the numbers in Table | are given for A= 1. 

Our analytical expressions for the form factors contain some differences with respect to Ref. as explained in 
Appendix |b| ||i7[] . Numerically they end up being most relevant for i — > 07. Fig. y illustrates the case in which a 
common value is taken for all the scalar masses, as might be useful for comparison |l8| with Fig. 2 of Ref. §. We 
can see that the analytical difference between us and Ref. Q translates into a numerical difference of more than one 
order of magnitude for the t — > C7 decay rate. 

From Table |, we see that Br{t 07), Br{t eZ), and Br{t — > eg) can be substantially enhanced with respect 
both to the SM and to the 2IIDM's with no FCSC's (i.e. model I and model II). Depending on the size of the FC 
couplings, in model III we can gain even more than two of orders of magnitude in each branching ratio. This is likely 
to make a crucial difference at the next generation of lepton and hadron colliders where a large number of top quarks 
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will be produced. Therefore these machines will be sensitive to signals from non-standard top decays and should be 
able to put stringent bounds on the new interactions involved, the FC couplings of model III in our case. In view 
of these future possibilities, a careful study of the FC couplin gs of the model is mandatory and we will analyze the 
constraints that emerge from existing experiments in Sees. VI -VII. 



V. TOP-CHARM PRODUCTION AT /i+^i COLLIDERS 



Another interesting possibility to study top-charm production is offered by muon colliders |g| . Although very much 
in the notion stage at present, a /i^/i^ collider has been suggested jl9|-|2^] as a possible lepton collider. Muon colliders 
are especially interesting for two main reasons. They can allow a detailed study of the s-channel Higgs and also they 
may make it feasible to have high energy lepton colliders in the multi-TeV regime. Neither of these goals is attainable 
with an collider. 

If muon colliders are eventually shown to be a practical and desirable tool, most of the applications would be very 
similar to electron colliders. One additional advantage alluded to above, however, is that they may be able to produce 
Higgs bosons (7i ) in the s-channel in sufficient quantity to study their properties directly 23 - 2^,|^ . The crucial 
point is that in spite of the fact that the pTH coupling, being proportional to m^, is very small, if the muon 
collider is run on the Higgs resonance, -^/i = m>^, Higgs bosons may be produced at an appreciable rate p^ , |23| -p5|,p[ . 



_At 



mjii the cross section for producing 7i, uu, normalized to (Tq = a^iX^ fi 



e^e ), is given by 



RiV.) 
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where is the branching ratio of Ti 



and ae is the QED ffiie structure constant. If the Higgs is very narrow. 



the exact tuning to the resonance implied in Eq. (g3|) 
production will then be given by |23] 



may not in general be possible. The effective rate of Higgs 



Rin) = 



arctan 



m-j-iS 



R{n) 



where it is assumed that the energy of the beam has a finite spread described by S 

m^(l - 6) < s < m^(l + 6) 



(24) 



(25) 



and s is uniform about this range. 

In a recent paper H we have considered the simple but fascinating possibility that such a Higgs, Ti, has a flavor- 



changing Titc coupling, as is the case in model III or in any other 2HDM with FCNC's. The process fi'^ii 



tc + tc 



will then arise at the tree level as illustrated in Fig. IJ. It will give a signal which should be easy to identify, is likely 
to take place at an observable rate, and has a negligible SM background. Thus the properties of the important Htc 
coupling may be studied in detail. 




H c n 

FIG. 4. Tree level contributions to n'^ fi 




tc + tc in model III. 



For illustrative purposes we take Ti = ft." in model HI where a = Q (case 1) or tt/A (case 2). The main distinction 
between the two cases is that in case 2 the decays ZZ, WW are possible while in case 1 they are not (see 

Appendix ^ for the relevant Feynman rules). Thus case 1 is very similar to Ti, = A^. This will matter in computing 



the total width of the boson, i.e 
calculation. 

In general the FC coupling of /i" to tc can be written as 



the Br{h'^ tc), while it was completely irrelevant in the e~'"e 
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FIG. 5. Rtc for (5 = 0, 10 ^ and 10 ^ in case 1 (set of solid curves) and case 2 (set of dashed curves). We also plot RiTL) in 
case 1 (dot-dashed) and case 2 (dotted). 



where xl and XR are in general complex numbers and of order unity if Eq. applies. In particular we will consider 
the case in which Act = A ~ 0(1) (see Eq. (p^) and xl and XR are real. We treated the more general case in Ref. |^ 
to which we refer for further details. 
The decay rate to tc is thus 



V{n^ tc) 



32TrM^ 



m, 



H 



\XR? + \XL? 



(27) 



and, r(7i tc) = r(7i ct) at the tree level that we are considering for now. 

As we did in Sec. [II for the e+e^ case, also in the /i^/i^ case we can define the analogue of i?*'^ in Eq. (nj) to be 



(28) 



where 5^ and B^ are the branching ratio for H ^ tc and Ti. ^ ct respectively. We estimate Rtc in the following two 
cases: 



i) case 1: a = ; 

ii) case 2: a — it/A . 

Using the result in Eq. ( p7| ) and taking the expressions for the standard partial widths for 7i (e.g. T{TL tt), 
T{H — > 66), . . .) from the literature [|ll[, we obtain the following results. In case 1, if m-^ is below the tt threshold, 
Rtc is about 10^^ — 1 and in fact tc makes up a large branching ratio. Above the ti threshold Rtc drops. For case 2 
the branching ratio is smaller due to the WW and ZZ threshold at about the same mass as the tc threshold and so 
Rtc is around 10~^. All these results are illustrated in Fig. ^, where we plot i?(7i) and Rtc with 6 = 0, lO^'^ and 10^^ 
in case 1 and case 2 . 

For a specific example, let us take mn = 300 GeV, i.e. co ~ 1 pb- For a luminosity of lO'^'^cm^^s^^, a year of 10^ s 
(1/3 efficiency) and for S = 10~^, case 1 will produce about 5 x 10'^{tc + tc) events and case 2 will produce about 150 
events. Given the distinctive nature of the final state and the lack of a SM background, sufficient luminosity should 
allow the observation of such events. 

If such events are observed, the collider offers the additional interesting possibility of extracting the values of 

the Xl and XR couplings in Eq. (Eq) separately. What is measured initially at a collider is Rtc- One is required 

to know the total width of the Ti, and the energy spread of the beam in order to translate this into T{7i ~> tc). This 
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FIG. 6. Tree level contributions to each different F'^ — F^ mixing, for F — K, Bd, D, in model III. 

then allows the determination of |xlP + IX-Rp (see Eq. (|27|)). To get information separately on the two couplings we 
note that the total helicity of the top quark is 

H* = -Hf-^^^^|^ (29) 

from which one may therefore infer \xl \ and \xr\- Of course the helicity of the t cannot be observed directly. However, 
following the discussion of ||2^,^ one may obtain it from the decay distributions of the top quark. Unfortunately in 
the limit of small nic the helicity of the c-quark is conserved. Hence the relative phase of and Xi? may not be 
determined since the two couplings do not interfere. 



VI. CONSTRAINTS FROM F° MIXING PROCESSES 



From the previous analysis we see that both e~^e~ , jJi'^ n~ — > ct + tc and t — > 07, Z, g could be of some experimental 
relevance depending on the size of the FC couplings of model III. As is well known from the literature on FCNC's, the 
most dangerous constraints on tree level FC couplings come usually from — F° mixing processes {F = K,Bd,D) 
[^-^. In these references we can find the bounds imposed on some tree level FC couplings by different F'^—F'^ mixing 
processes. Due to the specific structure of the couplings of model III and to the new phenomenology at hand, we 
think that a more careful analysis is due, which takes into account both tree level and loop contributions. We have 
examined the AF = 2 mixing processes in detail and concluded that both the K*^ — and the — mixings are 
particularly effective in constraining some FC couplings, while the experimental determination of the D^—D'^ mixing 
is, for now, not good enough to compete with the other two mixings. However, due to the different flavor structure 
of the — mixing, it would be extremely important to have a good experimental determination in this case as 
well. We will address the problem more specifically later on in this section. 

In a model with FCNC's, F^ — F^ mixing processes can arise at the tree level. Therefore they are likely to be 
greatly enhanced with respect to the SM where they appear only at the one loop level. Due to the good agreement 
between the SM prediction and the experimental determination of the mass difference in the — and i?° — B^ 
systems, any tree level contribution from elementary FC couplings needs to be strongly suppressed. This was the 
original motivation for imposing on the 2HDM's a discrete symmetry |9| which could prevent tree level FCNC's from 
appearing (model I and model II). Our goal will be now to verify if, in model III, the hierarchy imposed on the 
fjj' couplings by the ansatz in Eq. ( [lO| ) is strong enough to make the tree level contribution to any F'^ — F^ mixing 
sufficiently small to be compatible with the experimental constraints. 

For any F'^ — F'-' mixing we have evaluated, both at the tree level and at the one loop level, the mass difference 
between the mass eigenstates of the system, given respectively by 

MkAMk ^ Re{K''\{sd)v-A{sd)v-A\K°) 
MbAMb. ^ \{B°d\ibd)v^A(bd)v-A\m)\ (30) 
MdAMd ^ \{D'^\{cu)v-a{cu)v-a\D'')\ 

where we use the notation {qq')v-A = 97'^(1 — l^)<l' ■ The tree level contributions for each different mixing are shown 
in Fig. ^, while a sample of one loop contributions are illustrated in Figs. |^ and ||, for the box and the penguin 
diagrams respectively. 

In our analysis we have made some general approximations which we want to discuss first. 



i) We observe that, as in the calculation of the Ztc and "ftc form factors in Sec. Ill, also in this case the value of 



the phase a plays a minor role numerically. Therefore, as we did in Sec. [II, we will set a = in the following 
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FIG. 7. Box diagrams which contribute at one loop to the K —K and to the Bj^- 
case is obtained by appropriately replacing the external and internal quark states. 



s,b 

B°i mixing, in model III. The D°-D° 



analysis. Once that is done we can focus our attention only on the contributions coming from = h'^, H^ = A^, 
and The possibility of having FCNC's induced by and will give rise to the tree level AF — 2 mixings 
of Fig. |. 

ii) In the evaluation of the one loop contributions to the /S.F — 2 mixings, we will take advantage of the fact that 
they are due to scalar bosons whose couplings to the quark fields are proportional to the masses of the quarks 
involved and, for the charged scalar, to some combinations of CKM matrix elements. Therefore, in each case we 
will consider only the dominant contribution, which quite often will correspond to the diagrams with a heavy 
quark loop. This procedure is clearly more approximate than the exact calculation one uses to perform in the 
SM ||2^, but a similar order of accuracy would not be necessary in our case. Nevertheless, as a check of our 
calculation, we have reproduced the SM result for each mixing and used it as a reference point to fix the right 
relative signs and normalizations. 

iii) The AF = 2 effective interactions generated by the new scalar fields are often more complicated compared to the 
SM results pof , because the scalar-fermion couplings involve more chiral structures. Therefore, the evaluation 
of the mass difference in the various — systems will involve the matrix elements of operators other than 
just the SM one 

Oi^^ = O^LL = (/7^(1 - 75)9)(/7m(1 - 75)9) . (31) 



In general, matrix elements of the following operators are involved 

Of = 

= (/759)(/759) 

O'a = (/V759)(/7m759) (32) 
O^LR = (/V(l - 75)<z)(/7m(1 + 75)9) 
OL = (/(l-75)<7)(/(l-75)<z) 
OL = (/(l-75)<z)(/(l + 75)g) 



for F — K, Bd, D; f = s,b,c and q — d,u depending on the F*^ meson that we consider. The matrix element of 
Oaf=2 is usually given as 

(F«|Oi^^|F°) = Bp{F'^\Oi'^^\F°)yjA (33) 
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where the ratio of the matrix element itself to its value in the Vacuum Insertion Approximation (VIA) is 
expressed by the S-parameter Bp. Extensive non-perturbative studies of Bk and Bb exist in the literature, 
especially from lattice calculations p7[ . For our purpose, we just want to evaluate each matrix element in 
the VIA and use a common i3-parameter (the one for Oy^^)- Clearly this is only an approximation. Such 
an approximation may be problematic in the SM, where one aims to get a very precise prediction, but it 
suffices for our qualitative discussion. In particular, we will assume Bk = 0.75, Bb = 1, and Bd — 1 [ p7| . 
Moreover, according to Ref. and to analogous calculations we performed for those operators that were not 
considered there, we will use the following expressions for the AF = 2 matrix elements of the operators (for 
a = 5, P, V, A, VLL, VLR, LL, LR) in the VIA 









-\Mr^ 




(F" 


0?|F°)... - 






{F° 


0^|F").,. = 


2 , ,rO,F , 




(F° 


05|F°)w. = 






{F° 


Ovll\F'^)via 




- 


{F° 






ML = 


(F" 




' 3^ 


MEn = 


{F° 


Olr\F°)via = 





6 
6 



3 ^ 



(34) 



for F = K, Bd, D. All the previous matrix elements have been expressed in terms of the matrix elements of the 
only two operators which do not vanish on the vacuum, i.e. 



(35) 
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where Mp and mf indicate respectively the mass of the meson and of the quark of flavor /. The recent evaluation 
of the pseudoscalar decay constants {fp) can be found in the hterature [^9| . In our calculation we have used 
the following set of values: — 0.160 GeV, fs = 0.175 GeV, and Jd = 0.200 GeV, where the first value comes 
from experiments and the last two are representatives of lattice calculations. 

iv) No QCD corrections have been taken into account for model III. From the SM case [^,^ we know that these 
corrections can be important when a precise comparison with the experimental data is needed ]3^ ]. However, 
in model III they would affect the evaluation of the constraints at a much higher degree of accuracy compared 
to the approximations that we are adopting; therefore we neglect them. 

We now proceed to the discussion of the important results. In Eq. ( pO| ) we basically parametrize our ignorance of 
the FC couplings of model III introducing the A^'^ mixing parameters. Therefore, the constraints we are going to 
impose will allow us to deduce the order of magnitude of some of the A^''^. Due to the fact that the analysis of any 
po_pO jj^ixing will involve both tree level and one loop contributions, we will have to deal with several couplings at 
the same time. In order to simplify the analysis we will first take all the A^'^ parameters to be equal. Depending on 
the result of this first approach to the problem, we will consider the possibility that different couplings are differently 
enhanced or suppressed. According to this logic, we have considered the following cases, corresponding to three 
possible assumptions on the FC couplings of Eq. (|lO| ) . 

Case 1 : Ay ~ A common to all the FC couplings. 

Case 2 : A^i, Xdj <C 1 for i, j = 1, 2, 3, i.e. negligible FC couplings for the first generation and no assumptions on the 
other FC couplings. 

Case 3 : as Case 2 but with the further assumption that 

Abb, Asb > 1 and A^, Act < 1 . (36) 



The results for each F"— F° mixing {F = K, Bd, D) in case 1, case 2, and case 3 are reported in Table ^ where we also 
give the corresponding experimental results and the SM predictions . Both for the SM and for each 

case of model III we also specify whether the dominant contribution is due to tree level or to one loop diagrams. The 
different relevance of the tree level and one loop contributions and how this imposes constraints on some particular 
couplings will be explained as we go along the discussion in this section. Moreover, for each different choice of the 
couplings, we have varied ruh, tua, and nic in the range between 100 GeV and 1 TeV. This is the reason for the range 
of values that are given in Table for each case and for each mixing in model III. In this section we will discuss 
in particular case 1 and case 2. The scenario described by case 3 results from the constraints imposed by i?b, P and 



Br{B Xg^) and will be discussed in some detail in Sec. VI. For convenience, a summary of its results is reported 
in Table || as well. 

In Case 1, i.e. when all the FC couplings are parametrized in terms of a unique A, the leading contribution 
comes from the tree level diagrams of Fig. H. The one loop contributions are always subleading all over the mass 
parameter space. For each F^ — F^ mixing, there are two possible tree level contributions, mediated by an and an 
A'^ neutral field respectively. We want to consider them separately because in our analysis we will vary rrih and niA 
independently. Moreover the two contributions differ by the chiral structure of the resulting four fermion effective 
interaction, of S'-type for the exchange and of P-type for the A*^ exchange (see their Feynman rules in Appendix 
^). Since Mp 3> Mg , the distinction between the two tree level contributions becomes important and we write them 
as follows 

h° ^ AMp" = 2 U'/f) ' -^Mi (37) 



-Mf 

for F = K, B , f = s, b, c, and q = d or u depending on the mixing that is of interest. In particular, since 
Mp ^ Mg , we observe that an interesting possibility will be to have a light /i" and a much heavier A". In this way 
the A° contribution would not be too large, while we would still have a light scalar field with FC interactions. 

In Fig. ^we illustrate the constraints imposed on the parameter A by the tree level mixings in the K^—K'^, B^—B^, 
and D'^ — D'^ case respectively. We recall that A is common to the three FC couplings which govern the tree level 
mixings 
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Experiment 





One 


SM 


Loop 


Model III 


Tree 


(Case 1) 


Level 


Model III 


One 


(Case 2) 


Loop 


Model III 


One 


(Case 3) 


Loop 



3.51 • lO"^'^ 
1.4-4.6- IQ-^'^ 



18 in-l7 



10"^'' - 10 



20 in-19 



10"^" - 10 



3.26 ■ 10"^^ 
10"^^ - 10"^^ 



14 in-l3 



10""* - 10 



17 in-16 



10""' - 10 



< 1.32 • IQ- 

IQ-i'' - 10" 



17 in-i8 



IQ-^' - 10 



14 in-l5 



IQ-'-* - 10 



TABLE II. Experimental values and theoretical predictions in the SM and in model III for AAIf (in GeV) for different F^ — F'^ 
mixings {F — K,Bd,D). Case 1, case 2, and case 3 correspond to three possible scenarios in which different assumptions are 
made on the FC couplings, as described in the text. Each range is obtained by varying the parameters of the model over a 
large region of the parameter space, compatible with phenomenology and with the assumption of weakly coupled scalar fields. 
The leading contributions (tree level or one loop) are also indicated in each case. 



V V V 

The three curves represent the upper bounds on A imposed by the present experimental results (see Table ||) for 
different values of the neutral pseudoscalar mass niA, when we fix nih — 200 GeV. Lighter values of nih would not 
give any relevant change to these bounds. Therefore, for each value of ttia, the region above a given curve is ruled 
out by the corresponding F'^ — F'^ mixing. As we can see, the most relevant role is played by K'^ — K'^ and B^ — B'^, 
which constrain A to be definitely smaller than unity, even for large values of niA (i.e. tua ^ 1 TeV). We have verified 
that if the experimental precision on D'^ — D^ were increased by one order of magnitude, this mixing would also start 
to play approximately the same role as K'^ — K'^ and B^ — B^, so that the three lines in Fig. ^ would then roughly 
collapse into one line. We thus see that the AF — 2 mixings put severe bounds on the magnitude of the ■Clj'^ couplings 
when we require all of them to be proportional to a common A. For A < 0.1 (as in Fig. ^ for tua ''^ 200 GeV) it 
becomes difficult for the FC couplings of model III to play a role in any process. For instance, the case of top-charm 



production that we discussed in Sec. HI as well as any rare top decay would be too suppressed to be of experimental 
relevance. On the other hand there is no good reason to assume that all the Xij of Eq. ( p^ ) are equal. 

Therefore, let us consider next the more general case in which the Ay parameters are different for each different 
coupling. This is indeed the scenario that we identified as Case 2. We do not have enough specific measurements 
to constrain all of them independently. However we can make some general remarks. We first observe that the tree 
level mixings constrain only three couplings to be necessarily small: ^f^, and In particular, we can take the 
point of view that the bounds on A in Fig. ^ do apply only to Asd, Xm, and Auc- In general, we can assume that 
all the FC couplings involving the first family, including also ^^j, are suppressed. In case 2, the tree level mixings 
are suppressed enough that some loop contributions might become important. We will consider a loop contribution 
to become important when it is clearly bigger than the corresponding SM prediction, because only in this case it 
becomes possible to deduce a strong bound on the new FC couplings. 

The one loop diagrams that are most likely to be relevant are those which involve charged scalars only, because 
the charged couplings can also contain terms that do not involve any of the suppressed couplings (see Eq. (g)) of the 
first-family. This is not the case for the penguin-like diagrams of Fig. ^, which indeed turn out to be negligible, except 
for very light and unrealistic values of the neutral scalar and pseudoscalar masses. We mainly have to focus on the 
box diagrams of Fig. |^ with charged scalars. Moreover, we notice that, if of comparable size, the most relevant box 
diagrams are going to be those with only one (rather than two) charged scalar boson because they are proportional 
to the product of two A,;j instead of four Ay. Thus, if quantitatively relevant with respect to the corresponding SM 
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contributions, they may be as effective as the tree level mixing in constraining the A's. 

For the B^—B^ mixing, the only sizable contribution comes from the diagrams with a top quark in the loop, because 
the scalar field couples to the quarks as given by Eq. (p^, i.e. proportionally to their masses. For K'^ — K'^ mixing, 
the diagram with a charm quark in the loop can be even more important than the one with a top quark, because of 
the CKM couplings involved. Finally, in the D^ — D^ case, the most relevant loop contributions come from the box 
diagram with a bottom quark in the loop. For very light rah and wla, the penguin diagrams generated by and , 
with a top loop, can be comparable. 

The results of this analysis can be read off Table |l[ In spite of the high degree of arbitrariness that we have in 
model III, we do not find any region of the parameter space in which the one loop contributions in case 2 are definitely 
larger than the SM prediction. In other words, once the first family is assumed to decouple, the AF = 2 processes for 
K°—K^, B^—B^ and D^—D^ place no further constraints on the remaining FC couplings (i.e. and Therefore 
this analysis tells us that case 2 is compatible with the existing experimental measurements of the Ai^ = 2 mixings 
as long as the second and third generation FC couplings of model III do not have a Ay much larger than 1. This is 
an interesting result as far as our predictions for <T(e+e^ — s- ic + ct) are concerned. Due to the constraints we cannot 
enhance this cross section in any substantial fashion, but at least it is not suppressed and the values shown in Fig. |^ 
are roughly correct, with some room for a small enhancement. Similar considerations apply to the case of the rare 
top decays that we discussed in Sec. IV. 

Other more complicated assumptions on the FC couplings of the second and third generation could also be in- 
vestigated. In order to consider only those possibilities that respect the best fit of the available phcnomcnologi- 
cal constraints, we will postpone the analysis of other scenarios after the discussion of additional constraints from 
Br{B Xgj), the p parameter, and Rb, in Sec. VII. As we have anticipated in presenting the results of Table ||, 
we will discuss in particular another scenario, the one we have denoted as Case 3, which could be of some physical 
relevance. 



VII. BR{B Xs-y), Ap, Rb AND THE CONSTRAINED PHYSICAL MODEL 



As it is the case for 2HDM's without FCNC's (model I and model II) [Q, so also in model III it is very difficult to 
reconcile the measured value of the inclusive branching fraction for B Xg^ |37 



Br{B Xs-f) = (2.32 ± 0.51 ± 0.29 ± 0.32) x 10" 



(39) 



with the experimental results for i?b, the ratio of the rate for Z ^ 66 to the rate for Z ^hadrons. The present 
situation 11 is such that R'^^^ > Rf^ (- 3cr) [|| 
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^expt ^ 0.2202 ± 0.0016 , Rf^ = 0.2156 (40) 

and the value of i?™^' seems to challenge many extensions of the SM pO|,^. However, several issues on the mea- 
surement of this observable are still unclear and require further scrutiny pip2| . It is plausible that the experimental 
situation will change in the future. Therefore we may want to consider both the case in which the constraint from 
R^^^^ is enforced and the case in which it is disregarded. A third crucial electroweak (EW) observable in this analysis 
is given by the p parameter, which turns out to be very sensitive to the choice of the mass parameters of any new 



physics beyond the SM. In a recent paper |41| we have studied the problem in detail, considering two major scenarios 
in which the constraint from Rb is either enforced or disregarded. In the following we will summarize the main results 
of Ref. and discuss them in the context of the more general picture of model III that we have been tracing till 
here. 

Let us first recall how the presence of an extended scalar sector and of the new FC couplings affects the theoretical 
prediction for B{B Xg^), Ap and i?f,. 



The SM result for the inclusive branching ratio Br{B X^^) |43 441 can be modified in order to include the 



contributions from the new scalar fields and we obtain that in model III 



Br{B -> X,7) T{b sj) _ 6a fiMR)f^ m2 _^ ir-W. .m2> 



^ = 77775 ITk ^ ^ h / {\cri^b)\' + 1^^ (m,)n (41) 

Br[B ~i Xceve) l^b^ceve) tymc/rab) V / 

where f{mc/mi,) is the phase-space factor for the semileptonic decay and F takes into account some 0{as) corrections 
to both B XcCVe and B —f Xsj decays (see refs. |4l| , ^ for further comments). Cj^\mi,) and c!}^\mb) are the 
Wilson coefficients of the two magnetic type operators that occur in model III, for arbitrary ^[^'^ couplings (see 



Ref. Ilj), 



Qr"'"^^ = ^rribsaf^'il ± i^)bF^, (42) 

OTT 



in the AS" = 1 effective hamiltonian evaluated at the scale p = mb- The approximate calculation of C^{nib) and 
Cj{mb) shows that in model III the Br(B Xgj) is always larger than the SM one This feature is very general 
and the enhancement of the Br{B — > Xs"f) depends on the assumptions we make on the new FC couplings and on 
the masses of the neutral and charged scalar fields. 

From the analysis of Ref. it is clear that the Br{B Xg^f) is very sensitive to any enhancement of the FC 
couplings. The neutral scalar and pseudoscalar contributions, involving truly new kind of diagrams with respect not 
only to the SM case but also to model I and model II 2IIDM's, are proportional to and The charged scalar 
contribution depends on the charged couplings and ^.^.j^. for m — c,t and dj — s,b Q. In particular, the really 



leading contribution arises from the diagram with a top quark in the loop and the relevant couplings will then be: 

(■u,r ^ 



^Ys^ and According to Eq. (^ they are expficitly given by 



= , a = &b+&b (43) 

^ts — £,ss^ts + £,sb^tb , £,tb = ^sb^ts + ^bb^tb ■ 

Therefore the final prediction for the Br{B Xsj) will depend on the whole set of FC couplings which involve 
the second and the third generation. A strong enhancement of any of the Xij would confiict with the experimental 
prediction in Eq. (|39|), unless some very specific assumptions on the other parameters (couplings and masses) are 
made Q. 

Let us now consider Rb, defined as 

_ rjz^bb) 

= r(Z^hadrons) ' ^^^^ 

Neglecting all finite quark mass effects (to, ~ 0) Q, the generic expression for T{Z qq), for q = b, c, . . ., can be 
written as 

Viz ^ qq) = ^^Mz [(A,,,.)2 + (A,,«)^] (45) 

where ae is the QED fine structure constant, Nc the number of colors, s„, the Weinberg angle and A^ ^.(/j) the chiral 
left and right couplings of the Zqq vertex. They can be written as the sum of a SM piece plus a correction induced, 
in our case, by the new FC scalar couplings of model III 
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\.Lm^^TLiR) + ^T^R) ■ (46) 

Since Rh is in the form of a ratio between two hadronic widths, most EW obhque and QCD corrections cancel, in the 
massless limit, between the numerator and the denominator. The remaining ones, are absorbed in the definition of 
the renormalized couplings a and Sy^ (cw), up to terms of higher order in the electroweak corrections |4S|j50| , |36|] . As a 
consequence, the Aq j^(^) coupHngs will be as in Eq. (^, with A^^^^^^ given by the tree level SM couphngs expressed 
in terms of the renormalized couplings a and Sw (cw). This feature makes the study of Rb particularly interesting, 
because the new FC contributions may be easily disentangled in the Zqq vertex corrections. Using Eqs. p5| ) and (|4^), 
we can express Ri, in terms of i?^" and i?^" as follows 

^ + + ^^^^ 

where 

for q = b,c [|T]. In Eqs. and (^8|), terms of 0{{A^™j^^)'^) have been neglected and the numerical analysis 
confirms the validity of this approximation. The vertex corrections A^* and A^]|™ in model III will depend on 
the new FC couplings and on the scalar masses. As explained in Ref. |4l], the dominant contributions to A^U^ are 
due to the charged scalar and tend to further decrease the SM result, pushing the theoretical prediction for Ri, in 
the wrong direction with respect to the current experimental result in Eq. ( ^0[ ) . The neutral scalar and pseudoscalar 
contributions would increase Rb, but they are very small. In order to enhance them we have to make very strong 
requirements on the couplings and masses of model III, which correspond to the scenario we identified as case 3. This 
scenario and its compatibility with the other constraints will be discussed later on in this section. 

Finally let us consider the p-parameter, i.e. the radiative corrections to the relation between Myy and Mz- We 
may expect that the W and Z propagators are modified by the presence of new scalar-gauge field couplings (see 
Appendix [A[) . In fact, the relation between M\y and Mz is modified by the presence of new physics and the deviation 
from the SM prediction is usually described by introducing the parameter po [ p4|j5^ , defined as 

Po = (49) 
pM| cos^ Ow 

where the p parameter absorbs all the SM corrections to the gauge boson self-energies. We recall that the most 



important SM correction at the one loop level is induced by the top quark [p6 52 



Within the SM with only one scalar SU(2) doublet Pq"^"" = 1. In the presence of new physics we have 

PO - 1 + Apr™ (51) 
where ApQ*^* can be written in terms of the new contributions to the W and Z self-energies as 

The determination of from FNAL allows us to distinguish between pg and p ~ 1 -|- ptop- From the recent 
global fits of the electroweak data, which include the input for mt from Ref. [p3| and the new experimental results on 
Rb, Po turns out to be very close to unity. For Rb—R^^^^ as in Eq. ( ^ ) and rrit = (174 ± 16) GeV, Ref. quotes 

Po = 1.0004 ± 0.0018 ± 0.0018 . (53) 

This result clearly imposes stringent limits on the parameters of any extended model. Using the general analyti- 
cal expressions in Ref. and adapting the discussion to model III (making use of the Feynman rules given in 
Appendix ^), we find that 
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ApQ™ ~ ^ [sin^a G(mc, to^, mn) + cos^a G{mc, uia, rrih)] (54) 
8v Stt^ 

where all the terms of order (Af^^/m^) have been neglected and we define 

r<( ^ 2 mlm\ ml mlmj^ ml rn\mjjj^ , rn\ 
G{mc,mA,mH,h) = 5 5- log ^ ^ — log ^5 \ 5 5 — log —5^ ■ (55) 

Therefore the choice of the set of scalar masses will be crucial in order to make Apo compatible with Eq. (p3|). 

The results of our analysis of model III |4l| ] indicate that there are two main available scenarios depending on our 
choice of enforcing _Rf, as additional constraint or not. As far as the assumptions on the FC couplings are concerned, 



they correspond to what in Sec. VI we called Case 2 and Case 3. However, new restrictions on the mass parameters 
have been imposed by the additional constraints we have discussed in this section. Therefore we will enlarge the 
definition of case 2 and case 3 to include also the bounds imposed on the mass parameters. The resulting two 
scenarios will constitute the two available physical solutions of model III. We will devote the rest of this section to 
illustrate them in detail and summarize their relevant features. 

1. If we enforce the constraint from R'^^^^^ (see Eq. (^0[)), then we can accommodate the present measurement 
of the Br{B Xg'^) (see Eq. (|39|)) and of the AF^2 mixings (see Table |l|) and at the same time satisfy the 
global fit result for the p parameter (see Eq. (^)) provided the following conditions are satisfied. 

i) The neutral scalar /i" and the pseudoscalar A° are very light, i.e. 

50 GeV < nih - < 70 GeV . (56) 

ii) The charged scalar 7J+ is heavier than and A" , but not too heavy to be in conflict with the constraints 
from the p parameter in Eq. (|5^). Thus 

150 GeV <mc< 200 GeV . (57) 

in ) The Cf^ couplings are enhanced with respect to the ones, as described by the pattern we identified as 
Case 3 

Abb > 1 and < 1 (58) 

Xsb > 1 and Act <C 1 . 



The choice of the phase a is not as crucial as the above conditions and therefore we do not make any assumption 
on it. Within this scenario i?b can be predicted to be less than 2a away of We refer the reader to Ref. ^] 

for more details. From the previous requirements on the parameter of the model, we understand that it is in 
general very difficult to accommodate the present value of R^^^^ in model III. However, if we assume that the 
FC couplings (namely, the A^- parameters) are arbitrary and dictated only by phenomenology, then it is still 
possible to find a very small region of the parameter space in which, in principle, model III is compatible with 
the important experimental results. The values of the neutral scalar and pseudoscalar masses are required to 
fit the narrow window of Eq. (^6|) and are very close to their experimental lower bounds. In order to increase 
them and still agree with R'^^^^ , we would need a heavier nic and this would be in conflict with Eq. (p7|). 

We recall that similar difficulties are present in model II as well [ p6[ |. However the important difference with 
respect to model II is that this analysis of the constraints of model III gives us some hints on the possible range 
of the new FC couplings. This can be used to explore interesting experimental consequences in FC transitions. 
Let us review some of the most important ones. 

As we can read in Table ||, also in Case 3 the main contribution to the mixing comes from the one loop 
diagrams. Contrary to the B^—B^ and K^~K^ case, we note that in the D^—D^ case the mixing can be much 
bigger than in the SM, although still a couple of orders of magnitude below the experimental bound. This is 
the only case in which, in a F'^ — F'^ mixing, a loop contribution of model III can be much bigger than the SM 
prediction. Therefore an improved experimental determination of the — mixing would be very effective 
in constraining the model parameters for case 3. 

Moreover the possibility of having large and ^£ couplings, as predicted by this scenario and allowed by the 
present constraint on B^ — B^, seems to be particularly interesting for the study of some rare B-decays (i.e. 
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B l+l- and B XA+l' 



extensively in Sec. VIII 



L0[ and of Z — s- 6s + hs. Due to its importance we will discuss this subject more 



Finally, surprisingly enough, we have verified that the cross section for top-charm production and the decay 
rate for the rare top decays that we discussed in Sec. are not suppressed even if the ^Yj couplings are. The 
contribution from the neutral scalars and pseudoscalar are clearly negligible and the final result is dominated 
by the charged scalar contribution. In this case, the analysis has to be extended with respect to the description 
we give in Appendix in order to include the complete expression for the charged couplings (see Eq. (^). In 
fact, the contribution from the charged scalar will be dominated by the Cihargod coupling of Eq. (^) instead of 
by the ^charged oncs as we assumed in Sec. [n iV and in Appendix The results of Appendix ^ arc still valid, 
provided some changes in the , £^^, ... couplings are made. 

2. If we disregard the constraint from B!^"^^ there is no need anymore to impose the bounds of Eqs. (|5^)-(|5^) 
and we can safely work in the scenario of Case 2, where only the first generation EC couplings are suppressed 

A„,Adj<l for i,j = 1,2,3 (59) 

in order to satisfy the experimental constraints on the F^—F^ mixings. We will assume the EC couplings of the 
second an third generations to be given by Eq. ([l0| ) with ~0(1). 
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FIG. 10. Br{B 

In this case model III predicts a Br{B 



Xg^) compatible with experiments at the 2o'-level 



GeV, as we can see in Eig. |10|. As soon as is not enhanced anymore, the contribution of t 



for mc > 600 
ic neutral scalars 

and pseudoscalar is completely negligible. Therefore, both the value of the mixing angle a and of the neutral 
scalar and pseudoscalar masses (to//, rn/i, and toa) are irrelevant. In particular, Eig. nO is obtained for a = 
and values for rrih and ijia resulting from the fit to Apo Pl| 



rriH-, mh <rnc < mA 



and iriA < rric < mn, ruh 



(60) 



We note that none of the previous scenarios would give an enhanced value of Bb, because in that case mA and 
mh would be required to be equal and light (see Eq. (|56|)). In fact, in this scenario model III predicts Bb to be 
slightly smaller than Bf^ (and Be slightly bigger than i?f^) Q. 

The main theoretical predictions for case 2 have already been discussed in many sec tions of this paper: top-charm 
production in Sec. H, rare top decays in Sec. IV and F" — mixings in Sec. VI, 
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Process 


SM 


2HDM's 
Model I,II 


Model III 
Case 3 


Experiment 








> 9- 10"^^ 


> 4.3- 10"^^ H 


BriB° ^ fi+fi') 

Br{B° ^ Xs^i+fI-) 
Br{Z 6s) 


~ 4- 10"^ 

~ 7 - 10"*' 
~6 - 10"* 


10-9-10-* 
10"^ -10"* 


(Amm ^ 1) 

10"^- 10-" 
(Amm ^ 1) 

10"* - IQ-*^ 


< 8.3- 10"^ S 

< 8.4- lO"*^ H 

< 2.5 ■ 10"^ Ji 
< 5- 10^^ |60f 

? 



TABLE III. Values of AMs, (GeV), Br{Bs m"^A'"), Br{Bs Xsfi'^fi'), and Br{Z bs) in the SM, in 2HDM's 
with natural flavor conservation (Model I and Model II) and in case 3 of Model III. Each range is obtained by varying the 
parameters of the corresponding model over a large region of the parameter space, compatible with phenomenology and with 
the assumption of weakly coupled scalar fields. The present experimental bounds are also given. The top mass is taken to be 
mt ~ 180 GeV. 



VIII. B°-B'i, Z ^bs AND SOME RARE B-DECAYS 



We have seen in Sec. VII that the present experimental measurement of i?b (see Eq. (^0|)) suggests a new pattern 
of FC couplings for model III, that we called Case 3, according to which the couplings would be enhanced with 
respect to the ones. In this section we want to examine the compatibility of this assumption with the — 
mixing and its phenomcnological implications for the rare decays: B° fi'^ fi , Bs /i^M Z bs (we 

denote in this way the sum of the final states bs and 6s). 

We have summarized our present theoretical and experimental knowledge of these processes in Table III , where the 
predictions of case 3 of model III are compared with the SM results and with the results of the 2HDM's with natural 
flavor conservation (model I and model II). The experimental situation is still uncertain for all of the processes under 
study and we report in Table III the existing experimental bounds. It is clear from Table HI that in these processes 



there are good chances that continued experimental search could show deviations from the SM predictions. 

As is the case for the other AF = 2 mixings that we have examined in Sec. VI, B^^—Bl is extremely important to 
constrain the corresponding FC coupling, i.e. 
bound |57| 



However, for the time being the experiments give us only a lower 



AA/r, > 4.3 X 10"^^ GeV or 



AM£ 



> 10 



(61) 



Specifying Eq. ( |37|) to the i?° case, from Eq, 
tells us that the^ 
of Sec. VII that this compatibility is realized in case 3 of model III 



|6l| ) we get a lower bound for the FC coupling. Therefore B^—Bg 
{, coupling need not be small, i.e. Xsb can be somewhat bigger than one. We know from the analysis 

Moreover, we have verified that the presence of 
an enhanced coupling does not represent a problem for non-leptonic decays of the b quark, in particular for those 
decays which arise at the tree level via b sec. In the spirit of this qualitative analysis, we ask that at the quark 
level the contribution from new physics to the rate T{b — > sec) is appreciably smaller than the corresponding SM one. 
This is indeed the case for Xgt < 40 , i.e. in a large range of values. 

It is remarkable that the B^—B^ mixing does not prevent from belonging to that small region of the parameter 
space of model III which is suggested by R'^^^^. Therefore we want to further investigate the phenomcnological 
consequences of case 3 of model III in some physical processes that more directly involve the coupling 



i w the rare 

B-decays B and i?" IJ-^I^^ and the decay Z bs. 

The phenomcnological relevance of the decay B^ ^ — > and in particular of Bd^s tias been pointed out in 

Ref. . Although the experimental measurement is still poor |5^ , |59[] , this is a rare but theoretically very clean B-decay, 
which is not affected by large QCD corrections. In model III it can arise at the tree level via the exchange of a neutral 
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FIG. 11. R''" as a function of ^/s for the two scenarios of case 2 (solid) and case 3 (dashed). 



scalar or pseudoscalar with FC interactions. We think that the possibility of having a large coupling prompts us 
to reconsider i?" A''''M~- As is the case for B^, the prediction for fJ''^f^~ in model III can be enhanced at least 

by a factor of 10^ with respect to the SM and to the 2HDM's with natural flavor conservation. However, the range 
reported in Table III is obtained for a i^^^ coupling given by Eq. ( p^ with A^^ ~ 1 and different enhancements of the 
coupling. As we said from the very beginning, we do not want to consider here the implementation of model III 
in the leptonic sector. Therefore we will take the number of Table III just as an indication of the possibility of new 
interesting signals from this rare decay. 

The B Xs^i^ case could be even more interesting. In fact, as we can see from Table [II, better experimental 
bounds [59|j60[ | exist and are nowadays only one order of magnitude away from the SM prediction, which is known to 
very high accuracy (including QCD corrections, long distance effects, etc.) ||6l| , |6^ . Therefore this decay could become 
a good constraint for the and couplings. In model III there are two possible contributions at the parton level: 
a one loop transition due to the one loop induced effective Zhs vertex and a tree level transition directly mediated by 
a non-standard neutral scalar or pseudoscalar (i?°, or The tree level contribution crucially depends on the 
order of magnitude of the coupling and on the assumption we make on the coupling. In particular, if is 
enhanced and m^, are not too heavy (as in case 3), this contribution could become very important for given 
by Eq. (10) with A^^~0(1). However the dependence on does not allow us to make strong statements. As we 



explained before, the one loop contribution in model HI does not depend on and can be comparable, for large 
and to the SM contribution. Indeed the effective Zbs vertices in the SM and in model III are given respectively 

by 



cw 47rs^ 



(62) 



9 



Cw 47rs 



w 



where ae is the QED fine structure constant, sw the Weinberg angle and we have used the notation j'^ fj — (l±7^)/2. 
Comparing the Ajf^ couplings we note that (A5f)sM — 0.1 while (Ajf^jMODg depend on and For instance, 
when Abb — 40 and Asb — 10, (A5f)MOD3 — 0.3, (A^)mod3 — 0.06 and the contribution of model III to the 
Br{B XsfJ-^ fi~') becomes comparable or even larger than the SM one. Therefore, this decay can play an important 
role in confirming the case 3 scenario of model III and the compatibility of the model with the present experimental 
prediction for i?b and the most important other EW constraints. 

Finally, the effective Zbs vertex also induces a non negligible rate for Z bs. This rate is given in the SM and in 
model HI respectively by 



T{Z bs) 



SM 



47rs2j, 



(Ai^) 



SM 



(63) 
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In the scenario of case 3, T{Z — s- fes)MOD3 can be almost an order of magnitude larger than T{Z — > 6s)sm- This 
is reported in terms of branching ratio in Table [II , where the previous rates have been normalized to the Z width 
{Vz — 2.49 ± 0.007 GeV). Much of the relevance of this decay mode in model III depends on the enhancement of 
the ^/^ couplings. Therefore, any experimental bound would be extremely effective. We illustrate in Fig. |ll|the cross 
section for the related process e~ sb ~\- sh normalized to the cross section for e^e~ — > 'y* ^ [i^ [i^ ^ i.e. 

. + y ^ (64) 

(7(6+6 ^ 7* ^ [T^ ) 

starting from values of the center of mass energy below the Z-peak up to ^fs — 1 TeV. The upper curve corresponds 
to case 3 and the lower curve to case 2. As we can see, there is a considerable enhancement in the scenario of case 3 
that we are considering in this section. Although the Z — s- 6s events are not as distinctive as the te-production 
events. Fig. ^ seems to suggest that the experimental situation may be favorable and more experimental effort in 
this direction appears very worthwhile. 



IX. OUTLOOK AND CONCLUSIONS 

In this paper and in Rcfs. |^,^,^ we examined the phenomenology of FCSC's that may occur in extended models. 
We strongly share the point of view with many that the extraordinary mass scale of the top quark should prompt us 
all to reexamine our theoretical prejudices against the existence of such currents, especially involving the top quark. 

A very simple extension of the SM with another Higgs doublet leads rather naturally to such scalar currents. 
The model has the nice feature that experimental information can be systematically catalogued and guidance for 
further effort can be sought. The model has important bearings for some key reactions: e'^e~(/i"^/x~) — > te + ie, cc; 
t c-^{g,Z): D"-D", and ^O-i?" oscillations; B{Bs) l+l^, B{Bs) l+l^Xs, and e+e'{Z) bs + bs. Continued 
experimental effort towards these can hardly be over emphasized. The model also has important implications for 
Z bb and we want to stress that it is extremely important to clarify the experimental situation regarding Z ~^ bb. 
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APPENDIX A: FEYNMAN RULES FOR MODEL III 



In this appendix we summarize the Feynman rules for model III, which are used in the calculations presented in 
the paper. We choose to work in the 't Hooft-Feynman gauge. 

Fermion-Scalar couplings 

We present the Feynman rules for the couplings of the scalar fields (neutral scalar), (neutral pseudoscalar), 
and (charged scalar), to up- type and down- type quarks, as can be derived from the Yukawa Lagrangian of 
model III (Eqs. (0)-(0)). Following the discussion of Sec. |l[ these are the Feynman rules we need in our calculation 
of Rb. 




<Ui 
^ [i^cKM-ci^(i +75) - e[;-^cKM(i - 75)] 

Although the Cj^'^ couplings are left complex in the above, in practice, in our calculation we assumed they are real, 
i.e. Cij'^ — ^dj^* ^ we were not concerned with any phase-dependent effects. 

Gauge boson-Scalar couplings 

Here is a list of the Z-boson, M^-boson, and 7 interactions with model III scalar fields. We report them in terms 
of scalar mass eigenstates, hP , A", and We always have to take note of the relations (see Eqs. (^ and 

(^) between the scalar mass eigenstates and (i?", , H^, H^) and use the fact that neither ZH'^H^ nor ZH^H^ 
couplings are present p^ , p^ . We note the absence at the tree level of vertices like A^Z'^Z'^ and A^W^Wi. 
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APPENDIX B: ONE LOOP FORM FACTORS FOR THE Ztc AND THE -ytc VERTICES 



In this appendix we give the complete analytical expressions of the one-loop form factors for the Ztc, jtc, and gtc 
((7 = gluon) vertices, defined by 



1677^ V mt mt J 



(Bl) 



where V = ^,Z,g and in the gluon case we understand that g = g°-T'^ for a = 1, ... 8, T" denoting the SU(3) color 

matrices. This parameterization is obtained by reducing the most general tensorial form of the A^^^ vertex using the 
Gordon's decomposition and the gauge invariance of the currents involved. 
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Referring to the definition of the mass scalar eigenstates in Eq. (||) , and using the Feynman rules of Appendix 
we have computed corrections both for a = and for aj^O . At the qualitative level of our analysis this makes a minor 
difference. Therefore, we decide to give in the following the analytical results in the case of q; = 0. In model III with 
a = 0, H'^, and coincide with the mass eigenstates H^, hP , and and all the new FC contributions come from 
the second doublet. Therefore the form factors A^'^'^'^\ B^'^'^'S\ . . . are calculated by summing up all the one loop 
corrections generated by the neutral scalar the neutral pseudoscalar A'-' and the charged scalar Following 
Ref. [0, we will write each single form factor as the sum of four different contributions, i.e. 

piv) ^ p(v) ^ p(v) ^ piv) ^ p(v) F = A,B,C,D (B2) 

where -F^^"* denotes the contribution coming from those diagrams in which an h!^ neutral scalar is exchanged, f''^^ 

and the same for the neutral pseudoscalar A^ and for the charged scalar and finally F^'' represents the 

mixed hP-A^ contribution (see Fig. |]). 

The choice of a = and of the previous notation should also help the comparison with an analogous calculation 
reported in Ref. Q, in which the authors computed the Br{t cZ) and the Br(t c-y). Indeed, in model III, the 
effective one loop couplings Ztc and jtc (calculated for ^ 0) enter in the same way both in the calculation of the 
rates Y(t — > cZ, 7) and in the cross section for e^e~ — ^ 7*, — s- ic + tc. We find that our results are often different 

from the ones reported in Ref. 0|. In particular, we confirm the results for -F^^"* and F^\ whereas we have different 
analytical results for F^'^ and F^^K 

As a proof of consistency of our results, we have checked that the sum of different sets of diagrams is divergence 
free and that the final one-loop vertices satisfy the right Ward identities. We recall that, for each different form 
factor i^(^) (for F = A, B,C, D), the pole cancellation is verified separately for F^^-* + F^^ + Fj^\ and F^K In 
this respect we note that if we use the results reported in Ref. |Q| , there is one case in which the cancellation of the 
divergent terms among different diagrams seems not to be verified, namely the term A^ + A^ + Afj in the A^ form 
factor ]6^ . 

We finally note that in the analytical computation of the Feynman diagrams, all the quarks except the top quark are 
taken to be massless. Non-zero light quark masses are kept only in the ^fj^ couplings, according to the assumption 
made in Eq. ([l0|). 

Our results for the form factors are described in the following, where, for the sake of comparison, we adopt the 
notation of Ref. 0|. The form factors are indicated as Fh, . . . (fo r F — A,B, C, D), where the upper index is dropped 



in order to slightly simplify the notation of Eqs. ( |Bl| ) and (B2). The distinction between the different vector boson 
(7, Z or g) is made at the level of the couplings. 
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where we have introduced the following definitions 

P''' = xmj -f (1 - x)ml 

7]^ = (x + xy + y'^)m'f + {1 — x — y)mf^ — xyq^ 



7 



h _ W + (1 - x)ml 
x'^ml -I- (1 — x)m\ ' 



(B3) 



(B4) 
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denoting by = s the mass of the gauge boson (physical mass or invariant mass depending on the process considered) 
involved in the top-charm production process. The couplings and are defined as the following two linear 
combinations of the original and couplings 



^'^--fiifct^^l:) and = 



(B5) 



while ay — cty^'^ and aA = ctji^'^ denote the vector and axial-vector couplings of the different gauge bosons, given 
respectively by 



-e 



(B6) 



4^ , \ [l-^sin^„, 
4 cos 9w \ 3 



4 cos 9v 







where e, gw and gs are the electromagnetic, weak and strong coupling constant, while 9w is the Weinberg angle. 
These results are in agreement with those of Ref . . The apparent difference in sign for Ch and Dh is only due to a 
different assumption on the momentum of the gauge boson. 
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(B7) 
(B8) 



where, in analogy with the previous case, we have introduced the following definitions 

= xm^ +{1- x)m\ 
rj^ = {x + xy + y'^)m1 + (1 — a; — y)m\ — xyq^ 

^ xmj + (1 — x)m\ 
^ x'^rnf + (1 — x)ra\ 

The couplings and are now defined to be 



(B9) 



and a-\&t{e.t+i. 



tc I 



(BIO) 



while av — cty^'^ and aA — oia^'^ as in Eq. (B6). In this set of results we find many points of difference with respect 
to Ref. §. 



Am = ctM / dxdy 
Jo Jo 
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log — 5- + m: 
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where we have used the foUowing definitions 

Ml /-I o I I 2\ 2 I 2 I 2 2 

r] = [1 — 2x — y + xy + X jm^ + yrrij^ + xm^^ — xyq 
?7 ^ ^ = (1 — 2a; — y + + a; jm^ + xm^ + ym^j — xyq 

The coupUngs are now expressed in terms of S,y\ and iy\, defined as 

fMi _^ fU (fU _ cU*\ tMi _ ^ fU (fU . f 
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(B13) 



while the couplings of the photon and of the Z boson to the neutral scalar and pseudoscalar are given respectively by 



2 cos 9v 



and = 



(B14) 



Also in this set of results we find many points of difference with respect to Ref. . 
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where the following definitions have been used 
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rp^ = (1 — a; — y)(mp — xni^) — xyq^ 
rf'^ = (a; + y)m% - y(l - a; - y)mf - 



a;yg 



(B15) 



(B16) 



Using Eq. (H), we should define ^'^ to be a linear combination of many terms. However, if we follow the Cheng and 
Sher ansatz in Eq. ( p^ with Ay ~0(1), there is only one term among all the possible ones which gives the leading 
contribution, and therefore we define 



^SttSct 



(B17) 



If different assumptions on the couplings of Eq. are made, this statement will need to be modified. However, in 
most cases, the main analytical results are still valid. 

Finally the remaining couplings are given for the different vector bosons respectively by 
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0'l=le , hl = -\e , JZ = e (B18) 
z 9wSnr ,z 9w I -, 2 2 \ jz 9w ,-. r, 2 \ 



■'W 



In our calculation we further assume that the couplings are real, as we are not concerned with any phase 

dependent effect, and symmetric for sake of simplicity. This amount to set in the previous expressions ^tc*=^ct- 
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